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I.  Introduction 


Protective  headgear  has  long  been  used  to  reduce  injury  in  such 
activities  as  construction,  mining,  fire  fighting,  motorcycling,  and 
sports  (including  football,  baseball,  and  hockey).  In  addition,  the  use 
of  protective  headgear  is  becoming  more  widespread  among  some  consumers 
who  were  traditionally  not  associated  with  helmet  wear  (e.g.,  bicycle 
riders) . 

Many  manufacturers  respond  to  this  demand.  Even  for  helmets  designed 
for  the  same  activity  there  are  large  variations  in  style,  materials,  com- 
fort and  cost.  There  are  at  least  10  manufacturers  who  offer  more  than  50 
different  models  of  football  helmets  alone.  The  prospective  purchaser  of 
this  equipment,  whether  football  coach,  fire  chief,  industrial  representa- 
tive, or  individual  consumer,  is  faced  with  the  monumental  task  of  choosing 
among  the  helmets  offered  for  a particular  activity.  To  further  compli- 
cate this  decision,  the  most  important  variable,  the  degree  of  protection 
offered,  is  the  least  well  understood  and  the  most  difficult  to  evaluate. 

There  exist  many  test  methods  for  evaluating  the  protective 
characteristics  of  the  headgear,  and  standards  have  been  proposed  to  main- 
tain minimum  performance  levels.  Some  of  these  are  promulgated  by 
authorized  government  agencies  (National  Highway  Traffic  Safety  Administra- 
tion, Occupational  Safety  and  Health  Administration,  Law  Enforcement 
Assistance  Administration);  some  are  promulgated  by  long  existing  voluntary 
standards  organizations  (American  National  Standards  Institute,  American 
Society  for  Testing  and  Materials);  others  are  promulgated  by  groups 
established  for  specific  purposes  (National  Operating  Committee  for 
Standards  in  Athletic  Equipment).  These  standards  are  presently  in  various 
stages  of  completion.  It  is  not  unusual  to  find  more  than  one  test  method 
for  the  same  type  of  protective  headgear.  Sometimes  helmets  rejected  by 
one  method  will  be  passed  by  another. 

These  difficulties  are  in  part  due  to  the  limited  information  which 
is  available  on  head  injury,  especially  internal  head  injury.  Another 
factor  is  that  the  relationship  between  a test  method  for  protective 
headgear  and  the  amelioration  of  head  injury  has  not  been  sati sfactori ly 
explored.  Also  the  importance  of  simulating  the  injury  environment  in  the 
test  method  is  not  well  established. 

The  purpose  of  this  report  is  to  summarize  the  difficulties  in  the 
evaluation  of  protective  headgear  and  to  suggest  the  existence  of  a 
relationship  between  test  method  and  injury  prevention.  First,  the  state 
of  the  art  concerning  the  mechanisms  of  head  injury  and  the  relationship 
of  these  injuries  to  the  modes  of  energy  input  are  briefly  reviewed.  Next 
the  general  guidelines  in  the  design  of  protective  headgear  are  discussed 
and  present  test  methods  for  helmets  are  summarized.  In  order  to  provide 
a quantitative  framework,  a simple  model  which  incorporates  many  features 
of  present  test  methods  is  defined  and  executed.  A rational  approach  to 
the  design  of  meaningful  test  methods  for  protective  headgear  is  duscussed, 
with  reference  to  the  results  of  the  model. 


Injuries 
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1 1 . Mechanisms  of  Injury 


A cross  section  schematic  of  the  head  and  its  contents  is  shown  in 
figure  1.  The  brain,  a viscoelastic,  inhomogenous,  and  anisotropic  material 
is  encased  in  the  skull  which  is  relatively  rigid  except  for  the  opening 
through  which  the  spinal  cord  protrudes,  the  foramen  magnum.  The  geometry 
of  the  skull  in  the  frontal  regions  is  marked  by  small  curvatures  and  bony 
protuberances  in  contrast  to  the  nearly  spherical  top  and  rear  (occipital) 
regions.  For  more  details,  the  reader  is  referred  to  references  1/  and  2/. 

The  types  of  injuries  that  are  likely  to  occur  when  the  head  is 
impacted  have  been  well  documented  ]_/ - 3/,  and  only  general  categories 
will  be  mentioned  here.  Among  the  ways  to  classify  head  injury  is  by  the 
type  of  material  damaged:  namely,  skull,  blood  vessels,  and  brain  matter. 
Skull  fracture  is  the  structural  failure  of  the  cranial  bone,  usually  due 
to  the  process  of  crack  initiation  and  growth  when  tensile  stresses  arise 
during  impact.  Linear  fractures,  for  which  the  resulting  skull  geometry 
is  unchanged  except  for  the  appearance  of  a crack,  are  the  most  common  4 /. 
Internal  bleeding,  when  it  accompanies  head  trauma,  may  be  found  either 
within  the  brain  tissue  (e.g.  petechial  hemorrhage)  or  between  the  brain 
and  skull  (e.g.  subdural  hematoma).  Nearly  all  such  injuries  are  potentially 
fatal  and  require  surgical  intervention  5 /.  The  injuries  to  brain  material 
are  usually  divided  into  contusions  and  lacerations.  Contusions  are  found 
on  the  outer  surface  of  the  brain  and  are  marked  by  a discoloration  of  the 
brain  material;  a laceration  is  a tearing  of  brain  material  often  found  in 
the  regions  of  the  bony  protuberances  of  the  skull  ]_/.  These  injuries 
may  occur  separately  or  may  accompany  one  another. 

When  the  impairment  of  function  is  immediate,  the  injury  is  usually 
diagnosed  as  a concussion.  There  are  about  as  many  definitions  of  con- 
cussion as  there  are  authors  on  the  subject.  Thomas  6/  reported  that 
concussion  is  - "characterized  by  immediate  and  transient  impairment  of 
neural  function,  such  as  alteration  of  consciousness,  disturbance  of 
vision,  equil  ibrium,  etc.,  due  to  mechanical  forces."  7/  He  restricts  his 
attention  to  cases  of  unconsciousness  and  suggests  that  the  ultimate  cause 
is  mechanical  input  to  the  brain  stem  axis  (by  stretching  or  shear  strain). 
However,  Ommaya  8/,  quoting  from  Walker  9/  argues  that  "concussion  involves 
not  only  the  brain  stem  but  other  loci  in  the  brain."  He  then  goes  on  to 
say,  quoting  now  from  Symonds  1_0/  : 

"Concussion  should  not  be  confined  to  cases  in  which  there  is 
immediate  loss  of  consciousness  with  rapid  and  complete  recovery, 
but  should  include  the  many  cases  in  which  the  initial  symptoms 
are  the  same  but  with  subsequent  long  continued  disturbances  of 
consciousness,  often  followed  by  residual  symptoms  - concussion  in 
the  above  sense  depends  upon  diffuse  injury  to  nerve  cells  and 
fibers  sustained  at  the  moment  of  the  accident.  The  effects  of 
this  injury  may  or  may  not  be  reversible." 

In  this  sense  the  alteration  of  consciousness  known  as  concussion  will 
include  those  cases  in  which  the  diffuse  disturbance  of  functional 
damage  is  accompanied  by  structural  failure  which  may  be  irreversible. 
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Relation  between  gross  head  motions 
and  tissue  level  injuries 

It  is  important  to  distinguish  between  the  local  disturbance  of 
function  or  structure  and  the  gross  head  motions  which  accompany  these 
local  disturbances.  These  local  disturbances  are  always  due  to  material 
deformation  which  may  be  quantified  by  field  variables  such  as  strain, 
stress,  or  pressure.  The  local  deformation  can  be  related  to  the  gross 
motion  of  the  skull.  Following  an  impact,  three  distinct  motions  occur* 

(figure  2) : 

1.  Skull  bending 

2.  Translational  component  of  the  acceleration  of  the  center 
of  gravity. 

3.  Rotational  acceleration  of  the  head  about  the  center  of 
gravity. 

The  latter  two  modes  of  motion  are  inertial  and  may  result  without  impact 
to  the  head  (say  in  the  case  of  whiplash).  The  skull  bending  mode  results 
from  the  contact  phenomena  between  a foreign  object  and  the  scalp,  trans- 
mitting a local  force  distribution  to  the  skull. 

Each  of  these  modes  give  rise  to  local  deformations  which,  depending 
upon  the  magnitude,  are  capable  of  producing  internal  head  injuries.  Only 
the  skull  bending  mode  contributes  to  skull  fracture.  Many  experimental 
and  analytical  studies  have  attempted  to  understand  the  mechanisms  by 
which  the  energy  of  the  gross  motions  is  transferred  to  the  local  tissue 
level.  In  the  rest  of  this  section,  some  of  these  results  will  be  briefly 
mentioned. 

While  all  modes  are  operative  in  most  impact  situations,  the  experimental 
and  mathematical  model  studies  tend  to  consider  the  rotational  motion  separately. 
In  many  animal  impact  studies  the  translational  component  of  acceleration 
has  been  singled  out  as  the  most  important  variable  for  measurement. 

Accordingly,  most  tolerance  criteria  for  head  injury  have-  been  referred  to 
the  translational  acceleration. 

The  mechanism  for  tissue  damage  by  the  direct  impact  mode  (skull  bending 
and  translational  acceleration)  has  been  qualitatively  explained  as  follows: 
when  the  skull  is  impacted,  pressure  gradients  arise  in  the  fluid-like  brain 
(as  in  any  accelerated  fluid-filled  container)  11/.  Positive  pressures 
develop  in  the  brain  at  the  site  of  the  blow,  and  negative  pressures  (normal 
tensile  stresses)  develop  at  the  opposite  side  of  the  brain  (contre  coup 
site)  12/.  Transient  pressure  waves  in  the  brain  and  stress  waves  in  the 
skull  also  tend  to  magnify  the  amplitude  of  the  negative  pressures  1 3/. 

It  was  postulated  that  the  large  negative  pressures,  which  appear  to  exist 
at  the  contre  coup  site,  may  cause  cavitation  in  blood  or  brain  tissue  11/-14/. 
However,  contracoup  injuries  are  often  observed  for  occipital  impacts  but 
rarely  for  frontal  impacts  15/,  a descrepancy  which  cavitation  theories  often 
fail  to  resolve  16/. 


*If  the  direction  of  the  blow  does  not  pass  through  the  neck  axis,  there  will 
also  be  a rotation  about  this  axis. 
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In  studies  with  water-filled  cadaver  skulls,  it  was  shown  that  the 
presence  of  the  foramen  magnum  enhanced  the  resulting  pressure  gradients  Y\J 
It  has  been  suggested  that  the  established  pressure  gradients  may  force 
brain  material  through  the  foramen  magnum,  causing  shear  stresses  to  develop 
in  the  region  of  the  brain  stem  6/.  Other  displacements  of  brain  material 
have  also  been  observed  to  occur  1J_/,  17/,  18/.  Locations  where  these 
relative  motions  are  large  are  thought  to  be  potential  sites  of  injury.  In 
particular,  Gurdijian  suggests  that  these  "relative  motions  seem  to  be 
the  most  likely  mechanism  for  frontal  and  temporal  lobe  contusions  rather 
than  negative  pressures"  1J_/. 

Figure  3 shows  the  Wayne  State  University  Concussion  Tolerance 
Curve,  first  published  in  1962  19/.  These  data  were  based  on  cadaver  skull 
fracture  in  the  short  duration  range  (because  of  the  clinical  association 
of  linear  skull  fracture  and  concussion)  and  on  human  volunteer  studies  for 
impacts  of  long  duration.  It  is  seen  that  the  average  translational 
acceleration  required  to  produce  the  injury  depends  upon  the  time  interval 
over  which  the  acceleration  is  applied.  Haut  et.  al  noted  that  such 
behavior  is  to  be  expected  for  viscoelastic  material  like  the  human  brain  20/ 
While  the  short-comings  of  this  tolerance  curve  have  been  rigoriously 
discussed  21/-  23/,  it  remains  the  most  widely  accepted  compendium  of  con- 
cussion tolerance. 

Several  indices  of  injury  potential  have  been  derived  from  the  Wayne 
State  data,  based  on  the  near  linearity  of  the  data  when  plotted  on  log-log 
paper.  Gadd  developed  the  severity  index, 


where  a is  the  acceleration  expressed  as  a multiple  of  the  gravitational 
acceleration  and  T is  the  duration  of  impact.  According  to  Hodgson  4/: 

"Gadd  derived  a critical  value  of  1000  (the  units  for  SI  are 
seconds)  as  being  the  threshold  of  danger  to  ^ife.  . . He  recog- 
nized that  the  Severity  Index  may  be  a function  of  contact  area.  . . 
Later  on,  Gadd  suggested  a critical  value  of  SI  equal  to  1500  . . . 
for  distributed  impact"  24/. 

a 

Arguing  that  the  SI  does  not  really  match  the  Wayne  State  data,  Versace  2 
developed  the  head  injury  criterion: 


where  t-,  and  to  are  those  times  between  0 and  T for  which  the  above 
expression  is  maximized.  Again,  "a  HIC  level  of  1000  is  assumed  to  be  a 
threat  to  1 ife"  4/. 


t2  • 2.5 
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Other  tolerance  criteria  have  been  based  upon  1 umped-parameter  models 
of  the  head  22/,  26/,  27/.  One  example,  the  maximum  strain  criterion  (MSC), 
is  shown  in  figure  4a  26 / where  parameters  were  chosen  to  match  the 
impedance  response  of  animal  or  cadaver  heads.  The  maximum  strain  in  the 
model  was  related  to  the  onset  of  discernable  structural  damage  in  experi- 
mental animals.  "The  tolerable  head  strain  for  humans  (found  to  be 
.006  in/in)  was  calculated  from  the  mathematical  model"  27/based  on 
scaling  factors  that  were  developed.  An  acceleration  tolerance  curve  was 
developed  and  compared  to  the  Wayne  State  curve  (figure  4b).  For  side 
impacts  a different  tolerance  curve  was  indicated  28/ . 

Ommaya  and  his  associates  1_6/,  29/ -34/  have  conducted  a large  part  of 
the  experimental  research  to  identify  the  contribution  to  injury  by  the 
rotational  component  of  acceleration.  This  mechanism  of  injury,  originally 
proposed  by  Hoi  bourn  35/,  involves  the  transfer  of  shear  stresses  to  the 
viscoelastic  brain  by  the  rotating  skull. 

In  whiplash  experiments  with  monkeys,  rotational  accelerations  were 
measured  and  tolerance  curves  for  concussion  were  derived  29/ , 30/  (see 
figure  5).  It  was  found  that  the  tolerance  levels  were  reduced  if  the 
resulting  head  rotation  was- produced  by  direct  impact,  and  suggested  that 
"skull  distortion  effects  appear  to  reduce  the  level  of  head  rotation 
required  for  the  production  of  cerebral  concussion  and  contusions"  1 6/ . 

When  animals  wearing  protective  helmets  were  impacted,  no  concussions  were 
observed  at  the  same  rotational  levels  16/.  Thus  it  was  thought  that  the 
gross  head  motion  of  skull  bending  contributed  to  the  shear  stress 
mechanism  of  injury  29/.  Reasoning  from  a scaling  law,  also  attributed  to 
Holbourn  (in  which  tolerable  rotation  levels  for  different  species  are 
inversely  proportional  to  the  brain  mass  raised  to  the  2/3  power),  the 
primate  data  was  extrapolated  to  humans  to  obtain  tolerances  of  rotational 
acceleration  for  long  duration  impacts: 

2 

7500  rad/sec  for  concussion  of  99%  of  population 

1800  rad/sec^  for  concussion  of  50%  of  population 

In  another  series  of  experiments  31_/,  32/  an  attempt  was  made  to 
isolate  the  two  inertial  gross  motions.  By  rigidly  securing  primate  heads 
in  pots,  the  heads  could  be  subjected  to  purely  .rotational  or  translational 
motions  (the  axis  of  rotation  was  not  the  center  of  gravity  of  the  head; 
therefore,  a component  of  translational  acceleration  was  present,  which  was 
comparable  to  that  of  the  translated  animals).  All  rotated  animals 
suffered  concussions,  while  none  of  the  translated  animals  were  concussed. 

On  autopsy  it  was  found  "that  a greater  number  of  lesions  occurred  in  a 
more  diffusely  widespread  symmetrical  manner  in  the  rotated  group,  whereas 
only  a few  asymmetrical ly  placed  focal  lesions  developed  in  the  translated 
group"  8 /.  This  does  not  suggest  that  such  focal  lesions  are  not  dangerous 
(all  are  potentially  fatal)  but  rather  that  they  occur  without  concussion 
(as  diagnosed  immediately)  in  the  purely  translated  group. 

To  further  quantify  concussion,  functional  damage  was  assessed  by  a 
somatosensory  evoked  response  (SER)  method  32/,  33/ : "all  animals  in  the 

rotated  group  exhibited  neurological  evidence  of  experimental  cerebral 
concussion;  in  contrast,  none  of  the  translated  group  showed  this  effect." 
At  the  local  level,  similar  changes  in  the  SER  of  individual  nerve  fibers 
were  observed  under  mechanical  loads  34/. 
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Mathematical  modeling  of  head  injury 

The  usefulness  of  mathematical  models  lies  in  their  ability  to 
quantitatively  relate  the  gross  head  motions  to  local  tissue  deformation.  The 
effect  of  the  various  input  modes  can  be  isolated  and  studied.  Of  course, 
they  are  limited  by  how  well  the  material  properties  of  the  skull  and  brain 
can  be  established.  At  this  time  the  mechanical  properties  of  the  skull 
(see  e.g.  36/,  37/)  are  more  reliably  known  than  those  of  the  structurally 
more  complicated  brain  (see  e.g.  38/,  39/) . 

If,  in  addition,  experiments  can  establish  tissue  level  failure 
criteria  (including  functional  failure,  say  by  experiments  similar  to 
those  of  Thibault  34/  described  earlier),  then  the  models  can  be  used  to 
establish  tolerance  limits  for  the  gross  head  motions.  This  procedure  would 
be  as  follows: 

1.  Perform  experiments  on  tissues  that  relate  the  severity  of 
damage  to  the  history  of  the  local  deformation.  This 
relationship,  expressed  in  functional  notation*  is 

ST  =^(D) 

where  D is  the  deformation  tensor. 

2.  Use  the  model  to  derive  another  functional  which  relates 
severity  to  the  history  of  the  gross  head  motion  X (itself  a 
function  of  position  and  time): 

SG^(X> 

Two  examples  of  these  functionals,  SI  and  HIC,  were  discussed  earlier. 

They  relate  the  severity  potential  to  an  integral  of  the  translational 
component  of  the  acceleration.  In  the  MSC  approach  26/  the  acceleration 
criteria  is  derived  from  a given  injury  index,  mean  strain.  However,  this 
relationship  is  between  two  gross  modes  (translational  acceleration  and 
skull  inbending)  instead  of  between  a local  injury  mechanism  and  a gross 
head  motion. 

The  mathematical  models  to  be  discussed  can  be  cateoorized  accordinq 
to  whether  or  not  skull  rotation  effects  are  included.  Some  of  the  studies, 
which  do  not  include  rotation,  have  modelled  the  brain  as  an  inviscid  fluid 
41/  - 47/  in  order  to  identify  areas  which  are  likely  to  sustain  large  pres- 
sure cFanges  (both  positive  or  negative).  They  are  also  useful  in  predicting 
sites  of  the  skull  at  which  fracture  is  likely.  These  models  range  from  the 


*A  functional  is  an  operation  which  maps  histories  of  one  variable  onto  a 
one  dimensional  space  (see,  e.g.  Truesdell  40/).  The  simplest  case  is 
where  only  the  present  value  of  the  variable  is  important,  in  which  case 
the  functional  reduces  to  a simple  function.  Therefore,  for  skull  fracture, 
where  failure  occurs  when  the  tensile  stress  exceeds  a certain  limit,  the 
local  severity  index  would  reduce  to  the  maximum  value  of  the  tensile 
stress. 
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rigid  skull  model  of  Guttinger  41/  to  the  sandwich  skull  model  of  Akkas  44/ 
and  include  the  effects  of  non-sphericity  45/  and  glancing  blows  46/.  The 
finite  element  model  of  Shugar  1_5/  incorporates  many  complicated  structural 
features,  including  the  foramen  magnum.  Goldsmith,  using  a physical 
model,  accounted  for  the  effect  of  the  neck  47/.  Because  of  wave  propaga- 
tion effects,  these  models  predict  large  values  of  negative  pressures  at 
or  near  (depending  upon  the  sophistication  of  the  model)  both  poles  (impact 
side  and  opposite  site).  However,  in  all  of  these  studies,  the  duration  of 
loading  is  very  small  (one  millisecond  or  less),  a situation  which  magni- 
fies the  importance  of  the  transient  fluctuations  48/ . Advani  and  Owings  49/ 
have  predicted  that  the  impact  duration  of  the  skull  against  a rigid  flat 
surface  is  at  least  3 msec. 

Other  investigators  have  accounted  for  the  low  shear  modulus  of 
the  brain  material  41/-  48/  in  order  to  identify  sites  of  large  shear 
strain  even  in  the  absence  of  rotation.  In  these  studies,  the  pressure 
distribution  results  were  similar  to  the  inviscid  studies  described  above. 
Owings  and  Advani  49/,  50/  with  an  elastic  brain  model,  predicted  large 
values  of  shear  strain  at  both  the  midbrain  area  and  the  brain  surface. 

When  viscous  damping  is  accounted  for  48/,  51/ , the  magnitude  of  the  shear 
strain  decreases  and  the  locations  shift  slightly.  Hickling  and  Wenner  48/ 
predicted  shear  strain  values  of  about  .025  rad  for  a 3 msec  pulse  with  an 
average  acceleration  of  about  150g.  This  is  on  the  same  order  as  the 
shear  strain  for  structural  failure,  .035  rad  38/ . 

In  the  above  studies,  longer  loading  durations  were  considered  (up  to 
20  msec).  As  the  duration  increases,  the  magnitude  of  both  the  negative 
pressure  and  the  shear  strain  decreases.  This  result  would  be  contradictory 
to  the  Wayne  State  Tolerance  Curve  (figure  3)  if  only  the  magnitudes  of 
these  quantities  were  responsible  for  brain  injury.  However,  the  results 
of  these  investigations  also  show  that  at  the  local  level,  there  is  also 
an  increase  in  the  time  duration  for  which  negative  pressures  and  shear 
strains  exist.  Therefore,  it  is  likely  that  the  injury  has  some  history 
dependence  at  the  tissue  level  itself. 

In  rotation  studies  the  skull  is  modelled  as  a rigid  sphere  which  is 
subjected  to  prescribed  rotational  accelerations,  either  step  functions  52/ , 
53/  or  finite  pulses  54/-  57/ . When  the  brain  is  modelled  as  an  elastic 
material,  Holbourn's  2/3  power  scaling  law,  discussed  earlier,  was 
confirmed  54/.  If  viscous  damping  effects  are  included,  the  magnitude  and 
transient  fluctuations  of  the  shear  strain  are  reduced  55/,  and  the  2/3 
power  scaling  law  applies  only  approximately  54/.  Bycroft  picked  a specific 
brain  location  for  concussion  injury  (non-dimensional  radius  of  0.3)  and 
calculated  the  value  of  the  maximum  strain  required  for  correlation  with 
experimental  whiplash  studies.  For  both  squirrel  and  rhesus  monkeys,  this 
strain  value  was  .05  rad.  Using  this  same  location  and  shear  strain  level 
for  human  concussion,  a threshold  level  for  concussion  was  calculated  to 
be  3500  rad/sec  . Liu,  Chandran,  and  von  Rosenberg  55/  obtained  similar 
results  using  a finite  difference  procedure.  Ljung's  model  56/  accounted 
for  the  effect  of  the  falx,  the  near  rigid  portion  which  separates  the  two 
hemispheres  of  the  brain.  Lowehhelm  57/ , using  the  results  of  this  model, 
calculated  a tolerance  curve  based  on  a local  injury  criterion  which  was 
strain  rate  dependent. 
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III.  Headgear 

The  preceeding  summary  of  injury  mechanisms  is  useful  for  suggesting 
performance  criteria  for  protective  headgear.  The  experimental  and 
analytical  studies  confirm  that  all  three  modes  of  the  gross  head  motion 
are  capable  of  producing  local  deformations  which  are  potentially  injurious. 

At  this  time,  no  single  mode  can  be  identified  as  the  most  harmful.  The 
tolerance  values  for  these  gross  motions,  where  they  exist,  are  quite 
sketchy,  but  some  indications  are  clear.  The  dependence  of  the  injury 
potential  on  the  history  of  the  gross  motions  (effect  of  input  duration) 
is  a factor,  especially  for  impacts  of  shorter  durations.  The  translational 
acceleration  indices,  SI  and  HIC,  account  for  this  effect  and  should  serve 
as  guidelines  until  the  state  of  the  art  improves.  Though  analogous 
rotational  indices  have  not  been  suggested,  the  necessity  of  reducing 
rotational  accelerations  is  obvious. 

One  approach  in  the  design  of  protective  headgear  is  to  attempt  to 
reduce  the  gross  head  motions  by  placing  deformable  material  between  the 
head  and  the  prospective  impact  surface.  Analogous  to  the  simple  spring- 
mass  model  of  Appendix  1,  this  is  expected  to  change  the  inertial  response 
even  if  no  energy  is  absorbed  by  the  helmet  material.  In  this  analogy, 
the  percentage  increase  in  the  duration  of  the  acceleration  response  is 
the  same  as  the  percentage  decrease  in  the  average  (and  maximum)  accelera- 
tion. According  to  the  injury  criteria  for  which  the  acceleration  is 
weighted  to  the  2.5  power,  this  changed  response  is  favorable.  Any  energy 
absorbing  features  of  the  helmet  system  would  serve  to  further  improve 
the  response  of  the  head.  Another  purpose  served  by  these  deformable 
materials  is  to  spread  out  the  applied  load  over  a large  area,  thus 
reducing  the  gross  head  motion  of  skull  bending*. 

Several  studies  have  listed  some  of  the  general  parameters  of  a helmet 
system  with  qualitative  descriptions  of  their  relationship  to  head  protec- 
tion 1_9 /,  58/-61  / . More  recently,  simple  one  dimensional  models  of  head/helmet 
systems  were  used  as  a guide  in  optimizing  helmet  design  parameters  62/-65/ . 

In  these  studies  the  helmets  were  modelled  as  Kelvin  elements,  a spring  and 
a dashpot  in  parallel.  In  Liu's  model  the  helmet  element  was  coupled  to 
a fluid  filled  cylinder  and  the  injury  criteria  was  "the  averaged  time 
spent,  at  contrecoup,  beyond  the  cavitation  pressure  of  the  fluid"  63/. 
McElhaney  et  al  65/  coupled  the  helmet  element  to  the  two-mass  head  model 
described  earlier  26/  and  used  the  maximum  strain  criterion  for  injury. 

This  paper  is  not  so  much  concerned  with  proper  design  aspects  of 
protective  headgear  but  rather  with  methods  of  evaluating  whether  or  not 
a helmet  offers  suitable  protection.  In  an  ideal  helmet  performance  test, 
one  would  subject  a "head"  to  a realistic  impact  load  and  measure  the 
response  of  the  head  to  determine  whether  this  response  is  within  tolerable 
limits.  The  problem  of  determining  which  responses  to  measure  has  already 


*As  there  are  few  reports  of  skull  fracture  among  helmet  wearers,  it  is 
assumed  that  most  helmets  sufficiently  reduce  the  skull  bending  mode.  There- 
fore, the  question  of  tolerances  for  skull  fracture  has  been  avoided  in 
this  report,  though  this  particular  injury  is  the  most  well  understood. 

For  a complete  list  of  references,  see  Mahajan  3/. 
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been  addressed.  While  it  would  be  best  to  measure  local  deformation 
quantities  which  are  related  to  injury  at  specific  sites,  the  present 
state  of  the  art  only  suggests  guidelines  for  tolerance  of  the  gross  head 
motions.  Nevertheless,  this  limitation  has  not  slowed  the  development  of 
helmet  test  methods. 

The  "head"  to  be  used  in  a performance  test  is  another  source  of 
concern.  Ideally  the  best  head-form  would  incorporate  enough  essential 
features  to  produce  a response  similar  to  the  human  head.  Efforts  in  the 
development  of  humanoid  head-forms  have  been  reported  66/ - 70/ . The 
resilient  components  of  these  headforms  were  adjusted  so  that  the  response 
would  fall  within  a range  of  static  and  dynamic  (linear  acceleration) 
responses  of  cadaver  heads  66/ , 67/ . Other  comparative  tools  such  as 
impedance  response  69/  and  frequency  analysis  70/  have  also  been  used. 
However,  there  remain  questions  as  to  the  reproducibil ity  and  durability 
of  these  headforms  71/.  In  most  test  methods,  therefore,  a metal  headform 
is  used. 

The  ability  to  simulate  a realistic  impact  situation  is  another 
desirable  feature  of  an  ideal  performance  test.  In  nearly  all  helmet  test 
methods,  a helmeted  headform  is  dropped  onto  an  impact  surface.  The  ration- 
ale for  this  method  is  based  on  two  observations  72/ : 1)  in  laboratory 

experiments, cadaver  head  responses  were  observed  which  were  nearly  indepen- 
dent of  the  body,  and  2)  the  duration  of  most  impacts  is  much  less  than  the 
neck  muscle  reaction  time.  Present  test  methods  incorporate  a wide  variety 
of  impact  surfaces  and  drop  energies  even  in  evaluating  helmets  for  the  same 
activity.  If  the  purpose  of  the  test  is  to  insure  that  potentially  injurious 
head  responses  are  avoided,  then  it  is  necessary  that  the  impact  conditions 
of  the  test  be  such  as  to  simulate  (or  at  least  be  related  to)  a real  life 
impact  situation.  In  many  activities,  one  impact  mode  predominates,  but  in 
some  (e.g.  football)  there  are  several  equally  important  modes. 

To  illustrate  the  variety  in  test  methods  for  protective  headgear,  a 
partial  list  of  present  methods  is  shown  in  table  1.  (See  references 
72/-  81/.)  Most  methods  attempt  to  account  for  the  time  duration  effect 
Tthat  the  head  can  withstand  higher  accelerations  for  short  times;  see 
Wayne  State  curve,  figure  3)  by  the  "dwell  time"  concept  (where  limits  are 
set  for  the  time  duration  for  which  the  acceleration  exceeds  prescribed 
values).  It  is  interesting  to  ask  how  these  criteria  compare  to  a bio- 
mechanical injury  index  like  SI.  As  an  example,  consider  the  criteria  of 
the  FMVSS  218  standard.  For  simple  acceleration  traces  (e.g.  curves  which 
approximate  haversines  or  triangles),  a curve  representing  a barely 
safe  acceleration  response  is  completely  defined  by  specifying  any  two  of 
three  rejection  criteria  (figure  6).  The  corresponding  SI  values,  also 
shown  in  figure  6,  are  all  extremely  high  and  depend  upon  which  rejection 
criteria  are  specified.  Of  course,  such  a comparison  may  not  be  relevant 
as  these  methods  often  do  not  claim  to  represent  realistic  impact  situa- 
tions. Rather  they  are  one  means  of  ranking  helmets  in  one  impact 
configuration . 

The  possibility  that  this  ranking  would  change  if  resilient  headforms 
and/or  impact  surfaces  were  substituted  for  their  metal  counterparts  is 
usually  not  addressed.  In  one  study, for  which  the  same  helmets  were 
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tested  by  two  different  methods,  there  was  little  or  no  correlation  in 
helmet  performance  as  measured  by  the  different  methods  82/.  In  addition 
to  these  differences  in  test  method  (headform,  impact  surface,  impact 
energy,  and  rejection  criteria),  other  effects  may  also  influence  the 
proper  ranking  of  helmets.  E.g.  would  the  helmets  which  perform  best  in 
a typical  drop  test  method  also  offer  the  best  protection  against  head 
rotations?  Does  the  size  of  the  impacting  object  make  a difference? 


IV.  Model 


Descri ption 

In  order  to  address  the  questions  raised  in  the  previous  section  and 
to  provide  a quantitative  framework  for  comparing  one  test  method  with 
another,  a simple  one  dimensional,  lumped  parameter,  model  is  utilized. 

The  headform,  helmet  and  impact  surface  are  all  identifiable  elements  and 
a broad  range  of  test  conditions  can  be  studied.  The  scope  of  this  work 
differs  from  previous  modelling  attempts  62/ , 67/ , discussed  earlier, 
in  that  our  purpose  is  not  to  identify  optimum  values  of  helmet  parameters 
for  preventing  a particular  injury.  The  simplistic  models  for  the  helmet 
and  head  and  the  limitations  in  the  available  tolerance  data  would 
severely  reduce  the  applicability  of  such  an  exercise.  Rather  we  are 
interested  in  comparing  headform  responses  (acceleration  vs.  time)  when 
test  conditions  are  changed.  The  model  is  shown  in  figure  7,  and  each  of 
the  elements  are  explained  in  turn. 

The  headform  is  modelled  as  either  a rigid  body,  to  represent  metal 
headforms,  or  as  the  two-mass  resilient  system  suggested  by  Stalnaker 
et  al  83/  (figure  4).  The  force  FH  between  the  masses  is  given  by  (see 
figure  7): 

Fh  = K [L  - (X1  - X2 )]  + C(XZ  - xp  (1) 

where  L is  the  initial  separation  of  the  two  masses.  The  parameters  of 
this  model  were  evaluated  to  match  the  driving  point  impedance  character- 
istics of  cadaver  heads  (groups  1 and  2 of  table  2)  26/.  In  these  studies 
the  mass  M]  was  subjected  to  prescribed  accelerations,  and  tolerance  curves 
were  developed  to  correspond  to  threshold  strain  levels.  However,  in 
impact  simulations,  the  mass  M]  makes  contact  with  some  foreign  object 
(possibly  the  inside  of  a helmet)  and  accelerations  are  recorded  for  the 
mass  M^  (either  by  placing  an  accelerometer  at  the  center  of  gravity,  for 
humanoid  headforms,  or  opposite  the  impact  site,  in  cadaver  experiments). 

In  order  to  determine  whether  these  parameter  values  give  reasonable 
impact  accelerations,  the  acceleration  of  mass  M2  in  the  model  was  compared 
to  the  response  of  cadaver  heads  in  an  experimental  impact  situation. 
McElhaney  et  al  27/  reported  the  resulting  accelerations  of  cadaver  heads 
which  were  impacted  by  a flat  rigid  object  of  known  mass  (10  kg)  and 
velocity  (11  m/s).  Humanoid  headforms  have  been  developed  to  match  the 
acceleration  responses  in  these  cadaver  experiments  67/.  When  the  head 
model  was  subjected  to  the  same  input  conditions  (equations  of  motion  and 
initial  conditions  to  follow),  the  maximum  acceleration,  Amax,  was  much 
greater  than  in  the  cadaver  experiments,  and  the  time,  tmax,  at  which  the 
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maximum  occurred,  was  much  sooner.  These  quantities,  amax  and  tmav,  could 
be  matched  when  the  parameters  in  group  3 (table  2)  were  used  in  trie 
model*.  Slattenschek  and  Tauffkirchin  23/  using  a one  degree  of  freedom 
model  also  chose  a similarly  low  value  for  the  spring  constant.  Also 
static  force-displacement  curves  of  laterally  loaded  cadaver  heads  (and  a 
representative  humanoid  headform)  suggested  that  the  value  of  an  effective 
spring  constant  would  lie  in  the  range  0.7  - 3.5  MN/m  67/.  Therefore, 
in  most  of  this  work,  the  humanoid  headform  will  be  represented  by  the 
parameter  values  in  group  3 of  table  2. 

Both  rigid  and  resilient  impact  surfaces  were  simulated.  In  the 
latter  case  the  deformation  of  the  surface  was  assumed  to  be  quasi-static**. 
Then  the  relation  between  the  force,  Fs,  and  the  displacement  can  be 
derived  from  the  contact  relations  between  elastic  bodies  85/ , 

Fs  = B(X3  - d)3/2  (2) 


(see  figure  8)  where  d is  the  initial  thickness  of- the  helmet  liner  material 
and 
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E,  v and  R are  the  elastic  modulus,  poisson's  ratio,  and  radius  of 
curvature  respectively  for  each  body  (subscript  1 for  impact  surface, 
subscript  2 for  helmet).  For  an  impact  surface  made  of  polyurethane, 

Ei  = 28  - 63  MN/m^  (4000  - 9000  psi)  86/.  This  is  much  smaller  than  the 
modulus  for  typical  helmet  shell  materials  like  polycarbonate, 

E2  = 6.9  X 10-3  MN/m^  (10°  psi)  86/.  For  a flat  pad,  R]  = =0  , and  helmet 
radius  of  order  10  cm,  a range  of  values  for  the  parameter  B can  be 
calculated: 


B — 10  - 25  MN/m 


3/2 


Melvin  and  Roberts  87/  report  the  results  of  experiments  with  typical 
helmet  liner  materials  which  were  loaded  at  high  strain  rates.  The  visco- 
elastic nature  of  these  materials  is  evident.  To  simulate  these  deformations 
we  will  consider  the  simplest  viscoelastic  lumped  parameter  elements:  the 


*The  masses,  Mi  and  M2,  were  arbitrarily  chosen  to  be  the  averages  of  the 
side  and  longitudinal  cases.  K and  C were  chosen  so  that  amax  and  tmax 
would  fall  within  the  range  of  cadaver  responses. 

**According  to  Goldsmith  84/,  such  approximations  are  valid  if  the  duration 
of  contact  between  the  impacting  objects  is  long  compared  to  the  time  that 
it  takes  for  elastic  waves  to  traverse  the  objects.  Wave  speeds  in  typical 
impact  surface  materials  like  polyurethane  are  sufficiently  large  to  allow 
this  approximation. 
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Maxwell  element,  a spring  and  dashpot  in  series;  and  the  Kelvin  element, 
a spring  and  dashpot  in  parallel.  The  equations  for  the  force  in  these 
elements  (derived  in  Appendix  2)  are  (refer  to  figure  8): 


Maxwel 1 : 


Kelvin: 


(3a) 

(3b) 


where  E and  n are  the  elastic  modulus  and  viscosity  for  the  helmet  liner 
material,  A is  the  area  over  which  the  load  is  distributed  (see  Appendix  2), 
d is  the  initial  thickness  of  the  helmet  liner  material,  and  a dot  above 
a symbol  indicates  differentiation  with  respect  to  time.  It  was  found 
that  the  Maxwell  model  is  more  capable  of  reproducing  the  load-deflection 
curves  that  were  reported  by  Melvin  and  Roberts  87/  (Appendix  2)  and 
therefore  this  model  will  be  used  more  extensively  in  this  work.  This  is 
in  contrast  with  other  previous  one  dimensional  studies  62/,  65/  in  which 
helmet  materials  were  modelled  by  Kelvin  elements.  The  range  of  values 
for  the  parameters  E and  n were  chosen  to  represent  a reasonable  range 
for  helmet  materials  based  on  the  work  of  Melvin  and  Roberts  87/  (see 
Appendix  2).  Therefore,  in  this  model,  the  characteristics  of  a single 
helmet  are  described  by  a pair  of  values,  (E,  n). 


The  governing  system  of  equations  is  completed  by  force  balance 
relations  for  each  of  the  masses  (figure  7): 

¥i  - fh  - F . <4> 


For  a large  rigid  impact  surface  (X3  = d,  X3  = 0)  the  system  (1),  (3a),  (4) 
and  (5)  represents  four  equations  for  four  unknowns  X-j , X2,  F,  F^.  The 
initial  conditions  are  X-|(0)  = 0,  Xq  ( 0 ) = V,  X^(0)  = -1,  MO)  = V,  F(0)  = 0, 
where  the  velocity  V is  related  to  the  drop  height  h by  V = /2gh*. 

If  a resilient  impact  surface  is  to  be  included,  F = Fs  (neglecting 
helmet  mass)  and  equation  (2)  becomes  the  fifth  equation  for  the  now  unknown 
position  X3.  For  the  case  of  a small  rigid  mass  M3  impacting  the  head 
another  equation  must  be  added 

«3*3  ■ F 

and  the  initial  conditions  on  velocity  must  be  changed  to 

^(0)  = X2(0)  = 0,  X3(0)  = -V 

Note  that  the  mass  of  the  helmet  is  neglected  in  this  system. 
Preliminary  experiments  indicate  that  the  response  of  the  headform  is 
nearly  independent  of  the  helmet  mass.  It  is  likely  that  this  is  the 
reason  that  present  test  methods  do  not  include  the  mass  of  the  helmet 
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when  calculating  the  input  energy.  The  rationale  for  neglecting  the 
helmet  mass  is  as  follows:  only  a small  fraction  of  the  helmet  mass  is 

contained  within  the  area  of  contact;  the  potential  energy  associated 
with  the  rest  of  the  helmet  mass  is  not  available  for  compressing  the 
liner  material  and  instead  is  dissipated  by  flexural  wave  motions  of 
the  helmet  material  outside  of  the  contact  reqion.  This  latter 
behavior  was  observed  in  high  speed  motion  pictures  88/. 

*'  These  equations  were  solved  numerically  on  a Univac  1108  computer 
using  a modified  Runge-Kutta  method  as  outlined  by  Gill  89/ . The  accuracy 
was  verified  by  comparision  to  simple  cases  for  which  the  exact  solution 
was  easily  determined  and  by  computing  the  same  results  when  time  step 
was  reduced  by  an  order  of  magnitude.  The  output  variable  of  most  interest 
is  the  acceleration  of  the  headform  (the  mass  M2  when  the  resilient  head- 
form  is  used  in  the  model)  expressed  as  a multiple  of  the  gravitational 
acceleration.  Once  the  acceleration  history  was  obtained,  the  injury 
indicies,  SI  and  HIC  (see  section  on  mechanisms)  were  computed  numerically 
by  fitting  the  acceleration  data  with  cubic  splines  90/.  The  accuracy  of 
this  latter  procedure  was  verified  by  comparing  the  results  to  that 
obtained  with  simple  acceleration  curves  for  which  SI  and  HIC  could  be 
calculated  analytically.  The  computation  runs  were  terminated  when  the 
acceleration  dropped  below  zero. 

In  the  proposed  model,  a test  method  is  considered  to  be  a 
combination  of  headform,  impact  surface,  drop  height,  and  a parameter 
for  evaluating  performance.  The  simulations  for  all  of  these  elements 
are  intentionally  simple.  It  is  assumed  that  enough  detail  is  included 
to  indicate  trends  when  comparing  one  test  method  to  another.  The 
results  are  expected  to  be  useful  in  guiding  future  experiments  and  are 
not  to  be  considered  conclusive  by  themselves. 

Results 


Variations  in  each  parameter  were  referred  to  the  following 
reference  system: 

Headform: 


Rigid: 
Res i 1 ient: 

Helmet: 

Maxwell  model: 


M]  = 4.54  kg  (10  lb) 

M1  = .23  kg  (0.5  lb) 

M2  = 4.31  kg  (9.5  lb) 

K =1.75  MN/m 
C = 7.7  X 102  N-sec/m 

E =20.7  MN/m2 
n = 2.07  X 10^  poise 
d = 2.54  cm 
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Impact  surface: 

Resilient:  B = 10  MN/m^ 

Drop  height:  h = 1.83  m (6  ft.) 

Each  parameter  was  varied  through  a representati ve  range  while  the  others 
were  held  fixed  at  the  above  values,  In  addition,  the  effect  of  a smaller 
impacting  mass  M3  was  also  computed.  These  results  are  shown  in  figures  8 , 
through  16.  In  order  to  examine  the  effect  of  different  test  methods  for 
a range  of  helmet  parameters,  some  limited  performance  measures,  amax,  SI, 
and  HIC  are  listed  in  tables  3 to  6. 

The  effect  of  variations  in  the  helmet  dashpot  parameter,  n , is  shown 
in  figure  8 for  several  headform/impact  surface  combinations.  It  is  seen 
that  the  effect  of  changing  n depends  upon  the  test  method  and  becomes  less 
important  as  more  resiliency  is  added  to  the  system.  Similarly  the 
effect  of  changing  test  method  is  less  important  for  the  smaller  values 
of  n.  Generally,  it  will  be  seen  that  the  effect  of  changes  in  any  one 
parameter  decreases  in  importance  as  resiliency  is  added  to  other  parts 
of  the  system.  These  indications  are  independent  of  the  parameters 
chosen  for  measuring  performance.  However,  it  is  seen  that  HIC  is  more 
sensitive  than  amax  to  changes  in  n (although  not  shown  in  the  figures, 

SI  is  about  10%  higher  than  HIC  and  behaves  similarly). 

Similar  results  are  shown  in  figure  9 for  variations  in  the  helmet 
spring  parameter,  E.  The  same  trends  apply  as  described  above  for  n. 

Note  that  changes  in  this  parameter  seem  to  have  a smaller  effect  on  the 
acceleration  response;  i.e.  percentage  changes  in  E are  less  noticible 
than  a similar  change  in  n.  Also,  for  changes  in  E,  HIC  is  less 
sensitive  than  amax 

In  figures  10,  11  and  12,  the  effects  of  resilient  headform  parameters, 
K,  C,  M-j  and  M2,  are  shown  for  both  rigid  and  resilient  impact  surfaces. 

As  expected,  the  effect  of  variations  in  each  parameter  is  more  noticible 
for  the  rigid  impact  surface.  Of  these  headform  parameters,  the  accelera- 
tion response  is  most  sensitive  to  changes  in  the  spring  constant,  K.  For 
values  of  K above  3.5  MN/m  the  response  is  not  very  different  from  that 
of  the  rigid  headform. 

The  effect  of  the  resilient  impact  surface  parameter,  B,  is  shown 
In  figure  13.  Again  variations  in  B are  more  noticeable  for  the  rigid 
headform,  and  the  effect  of  headform  changes  is  more  noticeable  for 
large  values  of  B. 

Figure  14  shows  the  effect  of  drop  height,  h,  for  several  headform/ 
impact  surface  combinations.  Again  the  aforementioned  trends  apply. 

Also  of  note  is  the  very  regular  dependence  of  both  amax  and  HIC  on  h. 

It  is  easily  seen  that,  regardless  of  test  method, 

amax(h)oC  hV2; 


1 '.at  is,  the  maximum  acceleration  (or  maximum  force)  is  proportional  to 
the  impact  velocity,  all  other  conditions  being  equal.  Such  simple 
relationships  were  not  noticed  for  the  other  parameters. 

In  figure  15,  the  effect  of  changing  the  mass  of  a rigid  headform  is 
shown.  Note  that  both  indicies,  ama„  and  HIC  decrease  with  increasing 
headform  mass.  Considering  the  simple  spring-mass  analogy  of  Appendix  1, 
this  result  should  not  be  surprising.  However,  the  compression  of  the 
helmet  liner  material  increases  with  headform  mass,  and  there  may  be  a 
greater  tendency  for  the  liner  material  to  bottom.  Bottoming  aside,  these 
results  may  indicate  a problem  for  those  test  methods  which  only  prescribe 
a value  for  the  total  impact  energy,  as  both  higher  headform  masses  and 
lower  drop  heights  tend  to  lower  the  acceleration  response. 

As  described  in  the  previous  section,  the  governing  equations  were 
modified  to  include  the  case  of  a finite-sized  missile  striking  a 
helmeted  head  (or  headform).  Figure  16  shows  the  effect  of  varying  the 
mass  of  this  missile  while  its  velocity  remains  fixed.  The  response  for 
M3  = «=  is  the  same  as  when  the  head  impacts  an  infinite  rigid  mass. 

Note  that  as  M3  decreases,  the  durations  of  contact  are  reduced  as  well  as 
the  maximum  accelerations. 

In  table  3,  the  effect  of  headform  and  impact  surface  are  shown  for  a 
range  of  helmet  parameters.  Here  it  is  seen  that  the  general  trends 
referred  to  earlier  (figures  8 and  9)  are  consistent  for  the  entire  range 
of  helmet  parameters  considered.  If  the  maximum  acceleration  is  used  as 
a performance  measure,  there  are  several  examples  where  the  relative 
performance  (of  one  "helmet"  compared  to  another)  reverses  order  when  the 
test  method  changes.  For  example  consider  the  helmets  represented  by  the 
following  pairs  of  parameters:  E = 13.8  MN/m^,  n=  3.45  X 10^  poise  and 

E = 20.7  MN/m^,  n = 2.76  X 10^  poise.  The  former  performs  better  when 
tested  with  a rigid  headform  and  rigid  impact  surface;  as  the  resiliency 
of  either  test  method  parameter  increases,  there  is  a point  at  which  the 
latter  helmet  performs  better.  These  reversals  occur  because  the  quantita- 
tive effect  of  changes  in  test  method  is  not  independent  of  the  helmet 
parameters.  The  occurance  of  such  reversals  always  involves  a change  in 
both  of  the  helmet  parameters  (E  and  n).  For  the  examples  where  these 
reversals  occur,  the  computed  accelerations  are  not  very  different;  however, 
it  must  be  emphasized  that  only  a simple  helmet  model  is  being  utilized 
and  that  an  increased  likelihood  of  such  reversals  should  be  expected  if 
real  helmets  were  tested.  Experimental  evidence  of  such  reversals  in  the 
order  of  helmet  performance  have  also  been  reported  82/.  Note  also  that 
the  occurence  of  these  reversals  appears  to  be  far  less  likely  if  the  HIC 
or  SI  is  used  as  a measure  of  performance  rather  than  the  maximum  accelera- 
tion. In  fact  there  is  only  one  example  of  a HIC  reversal  in  table  3. 
Probably,  this  is  because  these  indicies  tend  to  spread  out  the  performance 
data  for  the  most  sensitive  helmet  parameter,  n . 

From  table  3,  it  appears  that  a rigid  headform/resil ient  surface 
combination  can  be  found  such  that  the  performance  measure  closely 
resembles  that  for  a resilient  headform  for  the  entire  range  of  helmet 
parameters.  That  is,  within  the  confines  of  this  model,  the  resiliency 
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of  the  impact  surface  can  be  used  to  compensate  for  the  resiliency  of 
the  headform.  As  seen  from  figures  10  and  13,  the  effect  of  changing 
the  impact  surface  parameter,  B,  is  similar  to  the  effect  of  changing 
the  most  important  headform  variable,  K.  Thus,  the  response  for  the  rigid 
headform/B  = 10  MN/m3/2  system  is  close  to  that  of  the  resilient 
headform/rigid  surface,  and  the  response  of  the  rigid  headform/B  = 5 MN/m3/2 
system  is  close  to  that  for  the  resilient  headform/B  = 10  MN/m3/2.  These 
indications  are  independent  of  the  performance  measure,  and  there  are  no 
examples  of  reversals  in  these  cases. 

Similar  results  are  presented  in  table  4 for  a range  of  spring  and 
dashpot  parameters  that  are  applicable  to  a Kelvin  model  (see  Appendix  2). 
This  model  gives  unreal istical ly  large  initial  accelerations  (because  the 
dashpot  is  in  parallel)  for  the  rigid  headform/rigid  impact  surface  and 
this  case  is  not  presented.  The  remaining  cases  show  the  same  trends  that 
were  evident  for  the  Maxwell  model:  1)  the  occurence  of  reversals  in  the 
order  of  helmet  performance  when  the  maximum  acceleration  is  used  as  a 
performance  measure,  2)  the  absence  of  reversals  when  the  HIC  is  used, 
and  3)  the  indicated  suitability  of  compensating  for  the  resiliency  of 
the  headform  by  the  resiliency  of  the  impact  surfaces.  It  is  also  of 
interest  to  note  the  close  agreement  between  the  two  models  when  the 
values  of  the  parameters  are  chosen  to  represent  the  same  helmet  liner 
material  (Maxwell:  E = 20.7  MN/m2,  n = 2.07  X 105  poise;  Kelvin: 

E = 2.07  MN/m2,  n = 1*38  X 105  poise;  see  appendix). 


The  effect  of  impact  by  a smaller  object  is  shown  in  table  5.  In 
this  table,  the  velocities  are  adjusted  so  that  the  impact  energy  in  each 
case  is  the  same  (in  figure  16,  the  velocities  were  constant),  and  the 
same  range  of  Maxwell  helmet  parameters  is  considered.  From  the  resilient 
headform  results,  it  is  seen  that  the  response  of  the  "head"  is  very 
sensitive  to  these  impact  conditions.  If  the  HIC  is  used  as  a performance 
measure,  the  impact  of  the  head  with  a large  object  is  much  more  severe 
than  if  the  head  is  struck  by  a small  object  even  though  the  energies  are 
the  same.  Again,  many  examples  of  reversals  are  evident  if  the  maximum 
acceleration  is  used  as  a performance  measure. 

To  further  illustrate  the  results  of  figure  15,  the  responses  for  a 
range  of  helmet  parameters  is  shown  for  two  cases  of  constant  energy 
impacts  in  table  6.  For  both  cases  a rigid  headform/rigid  surface 
system  is  used.  Large  differences  in  the  HIC  value  are  seen  for  the  two 
cases,  and  again  reversals  in  the  order  of  helmet  performance  are  indicated 
when  the  maximum  acceleration  is  used  as  a performance  measure. 


V.  Discussion 

As  the  preceeding  results  are  derived  from  a simple  one-dimensional 
model,  they  must  be  regarded  with  caution  until  confirmed  by  suitable 
experiments.  For  now  the  results  should  be  regarded  as  reasonable 
indications  of  expected  trends  when  test  method  parameters  are  changed. 
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The  real  life  impact  situation  should  be  the  starting  point  for  test 
method  development.  At  the  very  least  the  relative  velocity  at  which  the 
head  meets  an  impacting  object  and  the  properties  of  that  object  must  be 
described  for  one  or  more  modes  of  impact  which  pertain  to  a given 
activity.  A test  method  is  defined  by  a choice  of  headform,  drop  height, 
and  impact  surface  (input  parameters)  and  a measure  of  performance  (output 
parameter).  There  are  several  approaches: 

1.  Choose  input  parameters  to  simulate  the  real  life  impact 
situation.  The  tolerable  value  of  the  output  parameter  is  then 
based  upon  the  present  state  of  the  art  on  head  injury.  The 
model  shows  that  deviations  from  reality  in  the  headform  and 
impact  surface  are  capable  of  producing  reversals  (the  improper 
ranking  of  helmets) . 

2.  Compensate  for  changes  in  one  input  parameter  by  intentionally 
changing  another.  For  example,  the  model  suggests  that  the 
resiliency  of  the  headform  can  be  taken  info  account  by  an 
increased  resiliency  in  the  impact  surface.  In  such  cases, 
the  tolerable  value  of  the  output  parameter  would  remain  fixed 
at  the  appropriate  biomechanical  value. 

3.  Compensate  for  changes  in  an  input  parameter  by  adjusting  the 
output  parameter.  For  example,  if  a metal  headform  is  used, 
the  critical  value  of  the  performance  measure  could  be  relaxed. 

The  results  of  the  model  suggest  likely  trends  when  moving  from  the  real 
life  situation  to  the  test  method  situation. 

The  model  also  suggests  that  the  ranking  may  be  affected  by  the 
parameter  chosen  to  measure  helmet  performance.  The  present  state  of  the 
art  regarding  injury  mechanisms  strongly  suggests  that  the  likelihood  of 
injury  from  linear  accelerations  correlates  better  with  the  duration 
dependent  measures  like  HIC  and  SI  than  with  the  maximum  acceleration. 

The  results  of  the  parameter  studies  suggests  that  the  use  of  these  former 
parameters  is  also  less  likely  to  introduce  errors  in  the  ranking  of 
helmets  that  arise  from  changed  test  conditions. 

With  regard  to  evaluating  present  test  methods  the  model  can  be 
used  to  compare  present  rejection  criteria  to  the  biomechanical  parameters. 
As  an  example,  consider  the  rejection  criteria  of  the  FMVSS  218  motorcycle 
standard  (table  2).  From  table  3,  the  "helmets"  that  would  fail  for  the 
conditions  of  the  test  method  (rigid  headform/rigid  impact  surface),  are 
characterized  by  the  parameters:  E = 27.6  MN/m^;  n = 3.45  X 10^  poise; 

E = 20.7  MN/m2,  n = 3.45  X 10$  poise  and  E = 13.8  m/m2,  n = 3.45  X 105 
poise.  If  the  results  for  the  resilient  headform/rigid  impact  surface 
are  taken  as  representative  of  the  real  life  situation,  it  is  seen  that 
these  same  three  "helmets"  are  the  only  ones  for  which  the  SI  exceeds  1500 
(the  suggested  critical  value  for  distributed  impacts).  Therefore,  the 
practice  of  rejecting  helmets  when  the  acceleration  exceeds  a given  value 
for  a given  time  duration  may  be  an  acceptable  alternative  to  computing  the 
integral  measures,  provided  that  the  criteria  can  be  related  to  the  bio- 
mechanical measures  in  the  real  life  situation. 
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It  remains  to  discuss  some  of  the  limitations  of  this  work.  First, 
the  model  was  only  concerned  with  the  linear  acceleration  response.  If 
the  headform  were  allowed  to  rotate  in  some  realistic  manner,  the  likelihood 
of  a changed  helmet  ranking  is  a distinct  possibility  and  will  be  the 
subject  of  a future  paper.  Secondly,  all  of  the  conditions  considered  here 
were  based  on  an  impact  mode  in  which  the  head  is  relatively  free  to 
rebound  after  the  impact.  Therefore,  the  results  would  not  be  applicable 
for  a blow  to  the  top  of  the  head  (though  present  test  methods  do  not 
make  any  modifications  for  this  condition).  Finally  another  limitation 
is  that  the  helmet  model  has  no  provision  for  bottoming  (the  increased 
rigidity  of  the  helmet  liner  material  with  increasing  compression)  which 
is  sometimes  observed.  It  is  likely  that  the  effect  of  bottoming  would 
further  contribute  to  the  difficulties  in  properly  ranking  helmets. 

While  this  report  may  raise  more  questions  than  it  answers,  it  is  hoped 
that  a quantitative  framework  is  provided  which  can  be  used  as  a guide 
in  evaluating  how  well  a test  method  identifies  the  best  helmets  for  a 
given  activity. 
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Figure  4.  Head  model  (a)  and  tolerance 
curves  (b)  for  maximum  strain 
criterion  (from  Stalnaker, 
et.  al . 26/). 
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Figure  5.  Tolerance  curve  for  rotational 
accelerations  in  monkeys  (from 
Ommaya,  et.  al . 30/). 
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Figure  7.  Lumped-parameter  mathematical 
model  used  to  study  effect  of 
test  method  parameters. 
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Figure  10.  Effect  of  resilient  headform 
spring  parameter,  K: 

a)  rigid  impact  surface, 

b)  resilient  impact  surface. 
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Figure  11.  Effect  of  resilient  headform 
dashpot  parameter,  C: 

a)  rigid  impact  surface, 

b)  resilient  impact  surface. 
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Figure  12.  Effect  of  resilient  headform  mass 
parameters,  M-i,  Mo  (where 
M]  + M2  = 4.54  kgl:  a)  rigid 

impact  surface,  b)  resilient  impact 
surface. 
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Figure  13.  Effect  of  resilient  impact 
surface  parameter,  B: 
a)  rigid  headform,  b)  resil- 
ient headform. 
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Figure  15.  Effect  of  mass  of  rigid 
headform,  M:  a)  rigid 

impact  surface,  b)  resil- 
ient impact  surface. 
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Figure  16.  Effect  of  mass  of  a rigid 
missile,  M3,  striking  the 
helmeted  headform:  a)  rigid 

headform,  b)  resilient 
headform. 


Table  1.  Some  Features  of  Selected  Test  Methods  for  Protective  Headgear 
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Table  3.  Comparison  of  performance  measures  for  a range  of  Maxwell 
model  helmet  parameters.  E and  n . Entries  are  maximum 
acceleration  (amax),  severity  index  (SI),  and  head  injury 
criterion  (HIC). 
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Table  4.  Comparison  of  performance  measures  for  a range  of  Kelvin 
model  helmet  parameters  E and  n-  Entries  are  maximum 
acceleration  (amax),  severity  index  (SI),  and  head  injury 
criterion  (HIC) 
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Table  5.  Effect  of  mass  of  impacting  object  for  a range  of  Maxwell 
model  helmet  parameters  E and  n.  Entries  are  maximum 
acceleration  (amax),  severity  index  (SI),  and  head  injury 
criterion  (HIC).  The  case  for  M3  = ® is  the  same  as  for 
headform  striking  massive  surface. 
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Table  6.  Effect  of  rigid  headform  mass  for  constant  energy 
impacts  for  a range  of  Maxwell  model  helmet 
parameters,  E and  n.  Entries  are  maximum  accelera- 
tion (a  ),  severity  index  (SI),  and  head  injury 
criterion*  (HIC). 
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Appendix  1 


Consider  the  very  simplest  resilient  "head"  model  with  one  mass 
and  one  spring.  The  effect  of  a "helmet"  is  modelled  by  adding  a 
second  spring  in  series. 
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The  governing  equation  and  initial  conditions  for  this  system  are: 
MX  = -KX,  X(0)  = 0,  X(0)  = V 
The  acceleration  as  a function  of  time  is  given  by 


A sketch  of  the  acceleration  response  is  shown  below 


where  T is  the  time  at  which  the  acceleration  first  vanishes.  A severity 
index  can  be  computed  for  this  portion  of  the  acceleration  curve: 


Therefore  it  is  seen  that  the  severity  index  increases  with  K.  It  is 
also  of  note  that  both  the  maximum  acceleration  and  the  severity  index 
decrease  with  increasing  mass  of  the  "head." 


Appendix  2 


In  order  to  evaluate  the  proposed  models  for  helmet  liner  materials, 
and  to  determine  a realistic  range  of  parameter  values,  the  experimental 
results  of  Melvin  and  Roberts  87/  were  utilized.  In  this  study,  cylin- 
drical specimens  of  representative  helmet  liner  materials  were  loaded  at 
constant  velocities  and  the  compressive  stress  was  monitored.  The 
cross-sectional  dimensions  of  the  specimens  remained  constant  up  to 
maximum  strain.  A typical  stress-strain  curve  is  shown  below  (as  the 
solid  line  in  the  sketch): 


D 


CO 

CO 

CL) 

S- 

-l-> 

CO 

CL) 

> 

CO 

00 

<v 

i- 

CL 

E 

o 

o 


The  purpose  of  this  appendix  is  to  fit  two  simple  viscoelastic  models, 
the  Maxwell  element  and  the  Kelvin  element,  to  this  data  and  determine  the 
cooresponding  values  of  the  parameters  of  the  models. 

Maxwell  Kelvin 
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The  relationships  between  stress  and  strain  for  these  models  are  91/: 


Maxwel 1 : a + — a = Ee 

n 


Kelvin:  a = Ee  + ne 

For  loading  at  constant  strain  rate,  e = K,  the  Maxwell  equation  is  easily 
solved,  and  the  relationships  between  stress  and  strain  become: 

Maxwell:  c = nK  [l  - exp(-  e)] 


Kelvin:  a = Ee  + nK 


As  Melvin  and  Roberts  87/  noted,  their  experiments  were  conducted  at 
constant  velocity  rather  than  constant  strain  rate,  a difference  which  is 
unimportant  for  small  strains.  As  our  objective  is  to  determine  a 
reasonable  range  of  parameter  values  rather  than  an  accurate  curve  fitting, 
we  use  the  constant  strain  rate  assumption  up  to  maximum  strain  as  it 
simplifies  the  analysis. 

It  can  be  seen  from  the  above  equations  that  the  Maxwell  model  is 
much  more  capable  of  matching  the  stress-strain  curves.  The  Kelvin  model 
is  characterized  by  an  instantaneous  jump  in  the  stress  and  by  the  linear 
stress-strain  behavior.  As  shown  by  the  dashed  line  in  the  above  sketch, 
the  Kelvin  model  is  only  approximately  appropriate  if  the  early  non-linear 
portion  of  the  curve  is  restricted  to  sufficiently  small  strains. 

The  following  quantities  were  calculated  and  reported  in  the  Melvin 
and  Roberts  study  87/ : 1)  the  stress,  amax,  at  maximum  strain,  emax,  and 

2)  the  specific  energy  absorbed  up  to  maximum  strain,  defined  as  ■ 

xmax 

Fdx 

Q = 0 

Initial  Volume 

where  F and  x are  force  and  displacement.  If  A and  d are  the  initial 
cross  sectional  area  and  length  of  the  specimen,  Q can  also  be  expressed 
as: 
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For  the  two  viscoelastic  models,  the  expression  for  <^ax  and  Q are 


Maxwel 1 : 


°max  = nK  l1  * exp(-  nK  'rnax^l 
Q = nK{  E[1 


_ 1 

max  E 1 

L1  " exp(‘^ 

+ nK 

max 

Kelvin: 


max 


o — r max  [y 
Q - E — + nK  nmax 


Then  if  amax,  Q,  K and  emax  are  known  for  each  model,  there  are  two 
equations  for  the  unknown  values  of  the  material  parameters,  E and  n. 


The  data  which  corresponds  to  the  greatest  speed  in  the  Melvin  and 
Roberts  study,  217  in/sec,  is  used,  as  this  is  nearest  to  the  range  of 
impact  velocities  considered  in  this  report.  The  calculated  values  for 
polystyrene,  a typical  liner  material  for  motorcycle  helmets,  is  shown 
below. 
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As  both  E and  n,  for  both  models,  appear  to  vary  linearly  with  density*, 
the  values  for  more  typical  densities  (specific  gravities  of  about  0.1  for 
motorcycle  helmets)  can  be  extrapolated: 

Maxwell:  E = 20.7  MN/m2,  n = 2.07  X 10^  poise 

Kelvin:  E = 2.07  MN/m2,  n = 1-38  X 10^  poise 

The  range  of  values  to  be  considered  in  this  report  are: 

Maxwell:  E ~ 5 - 30  MN/m2 

n ~ 1 - 4 X 10^  poise 

Kelvin:  E ~ 1 - 4 MN/m2 

n ~ 0.5  - 1.5  X 10^  poi se 

To  derive  equations  (3a)  and  (3b)  of  the  text,  the  strain  e is  related 
to  the  displacements  X-j  and  X^  (see  figure  8)  by 

d - (X3  - X-,) 
e g 

and  the  force  is  given  by 
F 

0=  A 

f 

assuming  that  the  load  is  distributed  over  a constant  area**. 

*Such  linear  relationships  are  not  unexpected  for  polystyrene  foams  if  the 
volume  fraction  of  polystyrene  is  small  91/  as  it  is  for  the  densities 
considered  above. 

**In  a continuum  mgdel , the  contact  area  is  a function  of  the  force.  In  this 
one  dimensional  model  the  area  is  assumed  constant  and  equal  to  the  somewhat 
arbitrary  value  of  .016  m2,  chosen  to  give  realistic  headform  responses. 
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